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Bacterial colonisation of the gut is involved in the development of colitis-associated colorectal cancer. 
However, it remains unclear how the gut microbiota dynamically shifts correlating with colorectal 
carcinogenesis. Here, we reveal the longitudinal shifts in the microbial community that occur with 
colitis-associated colorectal cancer. High-throughput sequencing results for the bacterial 168 rRNA gene 
(V3 region) were compared for azoxymethane/dextran sodium sulphate-treated mice and control mice. We 
found that microbial community structure was significantly altered by chronic colitis. Microbes in the 
species Streptococcus luteciae, Lactobacillus hamster, Bacteroides uniformis and Bacteroides ovatus were 
increased during colorectal carcinogenesis. Histological measurements for a molecular network including 
six interconnected key factors from inflammation to cancer, namely p65, p53, COX-2, PPARy, CCR2 and 
P-catenin, indicated that the microbiome modifications were correlated with molecular pathogenesis of 
colitis-associated colorectal cancer. Phylotype modifications after each AOM/DSS cycle were identified. A 
longitudinal microbial network was then constructed for the gut microbiome and showed that the phylotype 
shifts during this process were complex and highly dynamic. This work may provide a deeper understanding 
of the role of the microbiota and microbe-host interactions in colitis-associated colorectal cancer. 

Experimental and clinical studies have frequently suggested that chronic inflammation is an important factor 
related to carcinogenesis'. Colorectal cancer is one cancer for which a relationship with chronic inflam- 
mation has been well established. It has been reported that individuals with inflammatory bowel disease 
(IBD) are at increased risk for colorectal cancer^'\ The mechanisms by which chronic intestinal inflammation 
promotes colorectal cancer have not been fully elucidated, but persistent exposure to inflammatory mediators and 
repeated damage/repair cycles are believed to contribute to the development of epithelial cell-associated cancer'''^. 

For cancers that are associated with inflammation, the microbiota is thought to be an important source of 
infection and a carcinogenic factor''. It was reported that some strains in Bacteroides fragilis and Escherichia coli, 
were enterotoxigenic and directly promoted tumourigenesis' ". In the colon, the microbiota is comprised of 10" 
prokaryotes and 10' different species'. The colonic microbiota is critical for mucosal tolerance under normal 
conditions and prevents pathogen infections'" ". However, recent studies have shown that the colonic microbiota 
can promote or aggravate chronic inflammation in the intestine". There is also increased evidence indicating that 
breakdown of the intestinal microbiota structure is associated with colon cancer'^ '"■. Thus, it is important to 
describe alterations to the intestinal microbiota during chronic intestinal inflammation and colitis-associated 
colorectal cancer. 

Structural changes to microbial communities are studied by various methods. One method that has been widely 
used in recent studies is the 16S rRNA gene sequencing approach. This approach provides an effective method 
with which to globally analyse a microbial community and identify specific phylotypes as biomarkers""'". Using 
the 16S rRNA sequencing approach, previous studies have shown that the gut microbiota structure is different 
between IBD patients and healthy controls'"". A more recent study that combined 16S rRNA sequencing and 
whole genome analysis of the colorectal cancer microbiota reported a modification in the gut microbiota structure 
in colorectal cancer patients and associated Fusobacterium with colorectal cancer^". These works provide large- 
scale descriptions of changes to the structure of the gut microbiota in IBD and colorectal cancer and demonstrate 
that gut microbial changes are closely related to these diseases, but further studies are still necessary to uncover the 
role the microbiota plays as gut inflammation progresses to colorectal cancer. 
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In addition to analysing clinical samples from IBD and colorectal 
cancer patients, animal studies have also been performed to uncover 
the pathological effects of gut microbial dysbiosis. Mice given 
chronic treatment of dextran sulphate sodium (DSS) treatment 
develop intestinal inflammation, and this treatment effectively mod- 
els IBD^'. By performing 16S rRNA gene sequencing and metatran- 
scriptomics surveys of DSS-treated mice, D. Berry et al. found a 
microbial community composition modification at the family level 
driven by DSS-induced colitis^^. Using a genetic deficiency mouse 
model together with the colon-specific carcinogen azoxymethane 
(AOM), J. C. Arthur et al. demonstrated that inflammation altered 
gut microbial composition, and Escherichia coli NClOl was specif- 
ically related to carcinogenesis". Recently, Zackular et al utilized a 
AOM/DSS model and suggested that the time course changes of gut 
microbiome were directly related to tumourigenesis^''. These studies 
show that animal models permit the effective study of microbial 
factors involved in the progression from IBD to cancer, and animal 
studies have revealed the important role of the gut microbiota in 
these diseases. However, details and a dynamic description of 
changes to the microbiota and the relationships among gut micro- 
biome, host molecules are still few as inflammation develops into 
cancer. The microbial link between gut inflammation and cancer is 
unclear. 

In this study, we employed a mouse model of colitis-associated 
cancer by combining DSS with AOM. By obtaining and sequencing 
mouse faecal samples after different cycles of AOM/DSS treatment, 
we developed a time-course study of the microbiota in colitis-assoc- 
iated cancer. First, histological measurements of p65, p53, COX-2, 
PPARy, CCR2 and P-catenin, six molecules closely related to inflam- 
mation and cancer were carried out. Based on the results, a dynamic 
molecular network was constructed to reveal the pathological pro- 
cess in the animal model of colitis-associated colorectal cancer. Next, 
gut microbial community structures were compared between control 
mice and mice treated with AOM/DSS. Differences in gut microbial 
composition were found between control and treated animals. Third, 
the phylotypes underlying these differences were identified. The 
species Streptococcus luteciae, Lactobacillus hamster, Bacteroides uni- 
formis and Bacteroides ovatus were significantly enriched in treat- 
ment group samples. The dynamic microbiome modifications were 
correlated with changes of tumour counts and the host dynamic 
molecular network. Finally, phylotype alterations were identified at 
different points along the time course. Longitudinal phylotype pro- 
files were mined to characterise bacterial fluctuations during the 
development of colitis-associated cancer. Some profiles contained 
microbes that may promote inflammation and cancer. A network 
of microbial profiles was constructed to further reveal the relation- 
ships between the dynamic patterns of phylotype variations. 

Results 

Progression from chronic colitis to colorectal cancer in a mouse 
model and a dynamic molecule network. BALB/c mice were treated 
with AOM and different cycles of 3% DSS. It was observed that the 
inflammation-induced tissue damage began after the first AOM/DSS 
cycle and became more evident after 2 AOM/DSS cycles (Fig. SIA 
and Fig. SIB). Number of dysplasia and carcinoma were counted in 
the colons of mice after treatment with three more cycles of AOM/ 
DSS (Fig. lA). Next, to characterise the progression of chronic 
colitis-associated carcinogenesis at the molecular level, immuno- 
histochemical staining for a molecule network consisted of p65, 
p53, COX-2, PPARy, P-catenin and CCR2 was carried out. These 
molecules were ranked high on inflammation and cancer related 
gene list inferred by CIPHER method^''. It also has been reported 
that these molecules are parts of main factors involved in colonic 
carcinogenesis^"'""', and they are also in close proximity to each other 
in a protein-protein interaction network (Fig. IC). In untreated mice, 
both p53 and p65 staining was almost undetectable. After AOM/DSS 



treatment, the expression level of p53 and p65 was increased and p65 
exhibited increased accumulation in the cell nucleus throughout the 
progression of carcinogenesis. The similar behaviors of p53 and the 
inflammatory factors were also detected in the DSS mouse model of 
colitis-induced cancer previously^'. P-catenin was increased and 
translocated from the cell membrane to the cytoplasm and 
nucleus. The expressions of COX-2 and CCR2 were increased 
during the development from inflammation to colorectal dysplasia 
and cancer. It is reported that COX-2 is highly expressed in the DSS 
model of colitis-related cancer^' and CCR2 is also activated in colon 
carcinomas^''. We also noted that the expression of PPARy in 
intestinal mucosa tissue and intestinal epithelial cells was elevated 
during this process (Fig. IB). And it was reported that in the absence 
of an exogenous ligand, PPARy in colonic epithelium has a 
protective effect against DSS-induced colitis, being embodied in 
the up-regulation of PPARy^**. The observed modifications of these 
molecules were quantified and visualized by dynamic molecule 
networks (Fig. IC). The dynamic patterns of these molecules may 
imply regulatory relationships among them. 

Diversity and structural changes of the microbiota community 
during colonic colitis-induced carcinogenesis. To study the 
diversity and structure of the intestinal microbiota community of 
the colitis-associated cancer murine model, 16S rRNA V3 region 
amplicon sequence libraries from faecal samples of AOM/DSS 
treated mice were sequenced. Sequencing gave approximately 22.9 
million paired raw reads in total. The number of raw reads pairs and 
the number of sequences represented OTUs of each sample were 
shown in Table SI. These high quality sequences were clustered 
into OTUs at different genetic similarity levels (100%, 97% and 95%). 

For alpha diversity analysis, rarefaction analysis was conducted for 
each sample. The rarefaction curves were constructed as the function 
of sequence counts to the number of OTUs observed (Fig. S2). 
Estimators of community richness, evenness and diversity were cal- 
culated based on the OTUs (97% similarity) (Fig. S3). We found that 
microbial richness was not changed by AOM/DSS treatment, while 
microbial diversity and evenness were significantly decreased in the 
model group mice with inflammation only (treated with one or two 
AOM/DSS cycles. Fig. S4). 

For beta diversity analysis, the Jaccard similarity coefficient was 
calculated based on the OTUs (97% similarity) to compare microbial 
community membership between control mice and all mice treated 
with AOM/DSS. Community membership was impacted by AOM/ 
DSS treatment, but the difference between the control and model 
groups was not significant (Fig. 2A, P-value = 0.07). The community 
membership of mice in the control group and mice in the model 
group with dysplasia and cancer was different, although this differ- 
ence was not significant after multiple test correction (raw P-value = 
0.046, P-value = 0.184 after Bonferroni correction). Next, to com- 
pare the microbial community structures of control mice and all mice 
treated with AOM/DSS, the Bray-Curtis similarity coefficient and 
the Yue & Clayton theta similarity coefficient were calculated. Both 
of these coefficients indicated that the microbial community struc- 
ture as significant changed by AOM/DSS treatment (Fig. 2A). The 
mice in the inflammation-only model group and the mice in the 
model group with dysplasia and cancer exhibited no significant dif- 
ferences in both microbial community membership and microbial 
community structure. Sample clustering under different distance 
metrics as visualised by non- metric multidimensional scaling 
(NMDS) methods (Fig. 2B). 

The relative abundances of phylum-level bacterial taxa were com- 
pared between control samples and samples from AOM/DSS-treated 
mice. As Fig. 3A and Fig. 3B shows, Bacteroidetes and Firmicutes 
were the dominant phyla in most mouse faecal samples, accounting 
for 22.7-85.2% and 8.99-62.8% of phyla, respectively, and 83.8- 
92.4% of phyla in total. Cyanobacteria, Proteobacteria, Tenericutes, 
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Figure 1 | Representative staining for molecular markers involved in the progression of chronic colitis-associated carcinogenesis. (A) Mice treated with 
1-4 AOM/DSS cycles were euthanized then the numbers of tumors in the colons were counted. Tumors were found in mice treated with 3-4 AOM/DSS 
cycles. Black bars are medians of each group. (B) Representative staining for six molecules reflects the pathological progression of chronic colitis and 
carcinogenesis. Rows represent the control group and mice treated with different AOM/DSS cycles. (C) The molecular markers examined are closely 
related in the protein-protein interaction (PPI) network. The dynamic changes in the expression levels (trans-nuclear rate for p65) of the molecular 
markers were quantified and represented by colour scale. Grey nodes are linker genes in the PPI and not for testing. 



Actinobacteria, Verrucomicrobia and Synergistetes for the most part 
constituted the rest of the gut microbiome in our samples (Fig. 3B 
and Fig. 3C). There were also some sequences that were bacterial but 
were not aligned to any taxonomic group by the methods we used. 
The less abundant phyla, i.e., with relative abundances < 10"''%, are 
shown in Fig. 3C and Fig. 3D. When comparing the relative abund- 
ance of phyla between the control group and the AOM/DSS-treated 
group, we found that the relative abundance of Bacteroidetes and 



Pirmicutes in the control group was not significantly different from 
the AOM/DSS-treated group (Fig. 3A). The other phyla that were 
consistently detected also did not change significantly. Among the 
rare phyla, Chlorobi was rarely present in the control group (6.25E- 
07% on average) but increased with increased numbers of AOM/ 
DSS cycles (3.44E-06% in mice treated with 1-2 AOM/DSS cycles, 
3.00E-05% in mice treated with 3-4 AOM/DSS cycles, shown in 
Fig. 3C). 
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Figure 2 | Shifts in microbial community structure induced by chronic AOM/DSS treatment. (A) AM OVA results for the between-sample Jaccard 
coefficient, the Yue & Clayton theta coefficient and the Bray-Curtis coefficient for the control group and AOM/DSS-treated group. The F-statistic is the 
ratio of the mean square among the groups and within the groups. (B) NMDS plots of samples with different distance metrics indicate the separation 
between control and treated groups. Mice with inflammation only and mice with tumours are not separated. Colours represent different mouse groups. 
Blue: control group. Red: Mouse group treated with 1-2 AOM/DSS cycles, only inflammation observed. Yellow: Mouse group treated with 3-4 AOM/DSS 
cycles, with observed dysplasia and carcinoma. 



Identification of phylotypes related to inflammation and cancer. 

As the above results show, the microbiota community structure of 
the colitis-associated cancer murine model changed, but phylum- 
level bacterial taxa were not greatly affected by AOM/DSS 
treatment. Therefore, an analysis was performed based on the 97% 
similarity OTU to determine whether a more specific phylotype 
difference existed between the control and the treated mice. After 
filtering out the OTUs with very low counts, 467 97% OTUs 
remained for further study. A multidimensional scaling plot 
(MDS) was drawn according to log-fold-changes distances between 
samples (Fig. 4A). The control group samples and AOM/DSS treated 
group samples separated clearly along the second dimension. Next, 
using an exact test, 70 97% OTUs were found to be significantly 
different between the control group and the AOM/DSS-treated 
group (Fig. 4b and Table S2, FDR < 0.05). Among these OTUs, 31 
were higher and 39 were lower in the model group treated with 
AOM/DSS. Two OTUs belonging to species Streptococcus luteciae 
(OTU01779, FDR = 0.015) and Lactobacillus hamster (OTU00907, 
FDR = 0.035) were significantly higher in the model group samples. 



Two OTUs frome genus Bacteroides, OTU00082 (classified as 
Bacteroides uniformis) and OTU00363 (classified as Bacteroides 
ovatus) were also significantly enriched in model group samples 
(Fig. 4B, Table S2); while three OTUs affiliated to genus 
Oribacterium (OTU00441, OTU00062, OTU00157) were enriched 
in samples from the control group (Table S2). These phylotypes 
could be considered as biomarkers to separate normal mice from 
mice with colitis or colon cancer. 

We also investigated whether there were phylotypes that were 
differentially enriched between mice with inflammation only and 
mice that developed dysplasia and cancer. The results showed that 
9 OTUs were more abundant in mice with inflammation only, and 
one OUT was more abundant in mice with cancer (Table. S3). Fig. 4C 
shows three of these OTUs that were most abundant in samples from 
model group mice. OTU00491 {uncliLSsiiied Erysipelotrichaceae) and 
OTU00024 (genus Bacteroides) were significantly abundant in the 
AOM/DSS-treated group (FDR < 0.05) and more abundant in mice 
with inflammation only than in mice with dysplasia and cancer. 
OTU00023 (unclassified S24-7) was significantly less abundant in 
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Figure 3 | Taxonomy classification of the mouse gut microbiome at the phylum level. (A) Relative abundances of the dominant phyla in samples of 
control group and model group. (B)-(D) Relative abundances of all phyla for mice in control group, group only with inflammation and group with 
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Figure 4 | Significant differential OTUs in the control group and AOM/DSS treated group. (A) Multidimensional scaling (MDS) plot, based on log- 
fold-changes distances between samples. (B) Log-fold-change and heat map of significantly differential OTUs (FDR < 0.05). (C) The three most 
abundant differential OTUs between mice with inflammation only and mice that developed dysplasia and cancer. *, FDR < 0.05, **, FDR < 0.01. "unci" 
means unclassified. In (A) and (B), each sample was denominated as follows: the beginning C/M represents the control or model group; the middle 
number represents the number of AOM/DSS cycles for the model group or the euthanisation batch for the control group; the last number represents 
different mice for repeats. 



SCIENTIFIC REPORTS | 4:4985 | DOI: 1 0. 1 038/srep04985 



6 



=3 

0 

c 

I 



E 
3 




B Cycle2 Cycles 

Cyclel 8X XIO ACycle4 




Cycle2 Cycles 
Cyclel ^ 2 X \ 4 V^Cycle4 



AOM/DSS cycles ' 




Figure 5 | The number of differential OTUs during the development of colitis-associated cancer. (A) The number of significantly increased OTUs (red) 
and decreased OTUs (blue) of mice treated with different AOM/DSS cycles compared to control mice (FDR < 0.05). Venn diagrams indicate the number 
of overlapping differential OTUs between mice with different AOM/DSS cycles: (B) increased OTUs, (C) decreased OTUs. 



mice with inflammation only than in control mice and mice with 
dysplasia and cancer (FDR = 1.71E-05) but exhibited no difference 
between the latter two groups. 

To identify which phylotypes were closely related to tumour 
occurrence, correlations between tumour counts and relative abun- 
dances of 97% OTUs (after filtering the low abundant OTUs) were 
calculated. 37 OTUs were observed significantly related to the 
colonic tumour counts of mice after different AOM/DSS cycles (P- 
value < 0.05, Table S6). Many of these OTUs classified as species 
Moryella indoligenes, Ruminococcus gnavus. Furthermore, the mole- 
cules measured in Fig. 1 reflected the developmental process of 
inflammation to colon cancer. We also explored correlations 
between changes of host molecular activity or expression with micro- 
bial community modifications. Some OTUs' changes were signifi- 
cantly correlated with host molecular expression or activity (Table 
S7). A dynamic network was constructed to visualize the coordinated 
changes of molecules and OTUs (Fig. S6). We observed that some 
OTUs significantly negatively related to the host molecule changes, 
such as OTUs affiliated to genus Oxobacter, species Roseburia faecis 
and Flexithrix dorotheae. These results suggest that some gut taxa 
were closely related to the development process of colitis to colorectal 
cancer and the dynamic patterns of these taxa could be important. 

Time course of differential OTUs across various AOM/DSS cycles 
and a dynamic network of OTU profiles. As Fig. 1 and previous 
studies indicated^'^'™, the development of colorectal cancer from 
inflammation is progressive and time-dependent. We next sought 
to reveal changes in OTUs after different AOM/DSS cycles across 
this progressive time course. After filtering OTUs with low 
abundance, a total of 19, 30, 48 and 73 OTUs were found to be 
significantly different between control mice and mice with AOM/ 
DSS treatment after 1, 2, 3 and 4 cycles (FDR < 0.05), respectively. 
Detailed information for these differential OTUs is shown in Table 



S4. The numbers of the increased and decreased OTUs are shown in 
Fig. 5A. Both the number of increased and decreased OTUs became 
higher as the number of AOM/DSS cycles increased with the 
exception of the number of increased OTUs after AOM/DSS cycle 
4, and this trend was more obvious for decreased OTUs. Further- 
more, both increased and decreased OTUs exhibited significant 
overlaps between successive AOM/DSS treatment cycles (FDR < 
0.05, Fig. 5B and Fig. 5C). This indicates that modification of the 
microbial community via AOM/DSS treatment is closely related to 
pathological progression. 

To stringently analyze the longitudinal microbiota community 
changes stringently, we employed STEM (Short Time Expression 
Miner)", a strategy which was previously used to analyze time course 
microarray expression data to mine the OTU profiles. First, a set of 
30 model profiles with five points was generated, and then each OTU 
was assigned to a profile. The 10 most enriched profiles comprised 
55.7% of all OTUs assigned to the 30 profiles. Three profiles were 
significantly enriched for OTUs that longitudinally varied during 
AOM/DSS cycles (Fig. 6A). Changes in OTU abundance for the four 
most enriched profiles (profile 28, profile 12, profile 20 and profile 
11) are depicted in Fig. 6B. Profile 28 is the most enriched profile; it 
exhibits a spike after 1-2 AOM/DSS cycles and a decrease to baseline 
or lower levels after 3-4 AOM/DSS cycles. More details about the 
OTU profiles are provided in Table S5. 

For most profiles, it was observed that OTUs in same profile were 
significantly associated in lineages than random background (Fig. S5 
and Table S8). We assumed that OTUs in the same profile shared 
overlapping niches. The OTUs in the same profile may represent a 
group of taxa which are functionally related and ecological coherent. 
We constructed a network to reveal the relationships between pro- 
files to explore the possible interaction between these bacterial 
groups (Fig. 7). The profiles showed different correlation relation- 
ships with each other. It was observed that some profiles with nega- 
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Figure 6 | Longitudinal profiles of OTUs reveal the dynamic patterns of microbial fluctuation during colitis-associated colorectal cancer. (A) Lines 
represent OTU relative abundance patterns across the five groups: the control group and the groups treated with 1-4 AOM/DSS cycles. The numbers in 
the upper-left corners are the serial numbers of the profiles. The numbers in the upper-right corners are the numbers of OTUs assigned to the profiles. 
The numbers in the lower-left corners are the raw P-values representing significance of overrepresentation. The numbers in the lower-right corners are the 
P-values adjusted by FDR. (B) The four profiles with the greatest overrepresentation; different OTUs are depicted with different colours. 
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Figure 7 | Microbial network reflects the relationships between the dynamic patterns of OTU profiles and the taxonomic constitution of profiles. 

Connections indicate positive (red) or negative (blue) correlations (Spearman correlation with P-value < 0.05). The size of each node is proportional to 
the number of OTUs in the profiles. Pie charts are overlaid on the nodes, representing the abundance of the family-level phylotypes in the profiles. 



lively correlated patterns demonstrated different phylotype consti- 
tutions, such as profile 20 vs. profile 10 and profile 26 vs. profile 12. 
This may indicate that bacteria within the same dynamic profile may 
be grouped together and interact with other groups of bacteria during 
colitis-induced colorectal carcinogenesis. 

Discussion 

The gut microbiome is a diverse and complex community whose 
balance is important for maintaining the healthy state of the host. 
Colonic dysbiosis is closely related to chronic intestinal inflam- 
mation and colorectal cancer'^. The transition from chronic inflam- 
mation to cancer is a multistep process"*^, and so the temporal 
association between the microbiota and disease development is key 
to understand the pathogenesis. However, there have been few stud- 
ies that have illustrated this temporal association. In this study, we 
employed the AOM/DSS murine model for colitis associated colo- 
rectal cancer and measured not only the overall shifts of the microbial 



community but also modifications to the microbiota over time and 
correlated with molecular changes. We found that microbial com- 
munity structure was changed by AOM/DSS treatment, and we iden- 
tified the phylotypes responsible for this change and associated with 
pathological process. When we generated a longitudinal view of gut 
microbiota modifications, it was found that phylotypes significantly 
changed after treatment with different AOM/DSS cycles, and 
dynamic patterns of the relative abundance of phylotypes were 
mined. 

AOM/DSS treatment exerts toxic effects on the colonic epithelium 
and alters the colon mucus barrier. AOM can induce DNA mutations 
and promote epithelial dysplasia. These effects result in intestinal 
niche alterations and expose the epithelial cells to pathogenic 
microbes. In accordance with previous reports that employed acute 
DSS treatment^^'", we found that the structure of the microbial com- 
munity was modified by repeated cycles of AOM/DSS treatment. The 
results showed that community structure was significantly changed. 
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but community membership was not (Fig. 2A). This indicated that 
the microbial divergence between AOM/DSS-treated and control 
mice was obvious in terms of relative abundance shifts of phylotypes 
but not for the absence or presence of phylotypes. We did not observe 
significant difference in microbial community structure between 
mice with inflammation only and mice with tumours. However, this 
result should be more carefully examined in future studies, as the 
sample size of these two groups in the current study was relatively 
small. 

At the phylum level, the taxonomic compositions of the gut micro- 
biota are similar between humans and mice. Bacteroidetes and 
Firmicutes comprise the largest portion of the entire bacterial com- 
munity'". We observed that the differences at phylum level were not 
significant. This result is in accordance with acute inflammation 
animal models^"'''"'. However, we found that the average relative 
abundances of these two phyla in our results {Bacteroidetes: 64%- 
69%, Firmicutes: 21%-26%) were different compared with prior 
studies {Bacteroidetes: 21%-31%^', 0.83%-9.3%"'; Firmicutes: 58%- 
65%^^ 86%-96%''*). These differences may be attributed to the ani- 
mal genetic background and experimental factors such as the loca- 
tion of sampling. The OTUs of the genus Bacteroides presented the 
same variation trends both in our work and in previous report^''. 
However, genus Odoribacter, which were increased after treatment 
in previous studies^^, exhibited a divergent change in our work. 
Genus Prevotella, which was decreased after treatment in previous 
report^"" showed opposite change in our work (Fig. 4B). Comparison 
of our results with previous reports suggests that the modification of 
some phylotypes with AOM/DSS treatment may depend on the gen- 
etic background of host animals, baseline structure of the microbiota 
and the AOM/DSS treatment protocols employed. 

Correlation analysis of microbial community modification with 
tumour counts and molecular expression/activity in the molecule 
network may imply the role of some taxa in the pathological process 
during colorectal carcinogenesis. For example, OTUs from members 
of species Ruminococcus gnavus were positively related to tumour 
counts, it implied a tumour promoting effect of the taxon. 
Accordingly, previous study showed that _R. gnavus could degrade 
secretory mucin and were increased in IBD"*^. While OTUOO 195 from 
species Roseburia /aec!s,which was negatively related to expressions 
of COX-2 and p53. R. faecis was reported to produce butyrate, a 
substrate important for maintain normal function of gut epithe- 
lium"". These results further suggested the intimate connection 
between bacteria and colitis-associated colorectal cancer; however 
more studies are stiU needed to explain underlying mechanisms of 
the molecule-microbe network for colorectal carcinogenesis. 

As the development of colitis-related cancer occurs in a stepwise 
manner, a time course study is helpful in understanding this pro- 
cess'^. Longitudinal OTU profiles were mined to describe the 
dynamic changes in the gut microbiota. The most enriched profile, 
profile 28, includes OTUs that were first increased and then 
decreased. These OTUs may be related to t inflammation process. 
For example, one of these OTUs, OTU00317, from genus 
Parabacteroides, strain from this genus was reported to have anti- 
inflammation effect"*. Another profile, profile 10, exhibited a spike 
after two AOM/DSS cycles just before carcinoma was observed 
(Fig. 6A). In this profile, OTU02959 (unci. Enterobacteriaceae, 
BLAST result showed it was from genus Escherichia, E-vlaue = 
4E-78). One strain in this genus, Escherichia coli NClOl, was 
reported to promote carcinoma in a mouse model of chronic colitis**. 
Another genus in profile 10, Streptococcus, was also found to be 
associated with colon cancer'''^. Additionally, we found that the genus 
Lactococcus was significantly enriched in mice treated with 3 and 4 
AOM/DSS cycles (Table S4). Both Lactococcus and Streptococcus are 
in the same family, Streptococcaceae, which was reported to be sig- 
nificantly enriched in human colorectal carcinoma samples^". 
Recently, a bacterial driver-passenger model was proposed for the 



involvement of the microbiota in colorectal cancer''". In this model, 
bacterial drivers were defined as the bacteria that induced DNA 
damage and initiated cancer. These bacteria may decrease in abund- 
ance after tumour formation, and some bacteria (called passengers in 
the model) that have a growth advantage can outcompete the drivers. 
Profile 10 in our result may represent a dynamic pattern of the 
bacteria drivers. 

There are complex interactions among bacteria in gut, and net- 
work models effectively descript such interactions'". In this work, we 
utilized dynamic profile clustering technique to explore the structure 
of microbial community. We found that the profiles were a group of 
taxa that were closely associated in lineages (Fig. S5 and Table SB). 
We assumed that taxa in such profiles may form coherent entities 
and have overlapping niches, just like co-occurring taxa''^. Taxa in 
different profiles could interact with taxa in other profiles as a whole, 
positively or negatively (Fig. 7). Thus, the microbial dysbiosis in 
colitis related colorectal cancer may be recognized as disorganization 
of interactions among different groups of taxa rather than individual 
taxa. But this hypothesis still needs more studies to prove. 

In summary, this work investigated the dynamic molecule- 
microbe networks related to colorectal carcinogenesis. By conduct- 
ing a longitudinal survey of gut microbiota modifications as well as 
molecular markers, we revealed a structural shift in the gut microbial 
community in a mouse model of chronic colitis-associated cancer. In 
addition to the finding that phylotypes that may play important roles 
in colitis and cancer and correlate with molecular pathological pro- 
cess, description of the phylotype profiles further provided a com- 
prehensive view of the dynamic changes to the gut microbiota that 
occur in colitis-associated cancer. This work may offer more 
information about the role of microbe-host interactions in cancer 
development. 

Methods 

Animal experiments. Male BALB/c mice between 6-8 weeks of age were used for 
experiments. Mice were kept in autoclaved cages and were given sterile food and 
water. The chronic colitis-to-colitis-associated colorectal cancer murine model was 
employed based on previous methods*''. The mice were divided into five groups that 
received different treatments. Mice in group 1 were given an intraperitoneal injection 
of AOM (10 mg/kg body weight). One week after the injection, the mice received 3% 
DSS (molecular weight: 36-50 kDa, Sigma, Germany) in their drinking water for 
seven days. The mice received another two weeks of sterile drinking water before 
euthanisation. Mice in group 2 received AOM and two cycles of DSS treatment as 
described above, mice in group 3 received AOM and three cycles of DSS treatment 
and mice in group 4 received AOM and four cycles of DSS treatment. Mice in the 
control group received sterile water before being euthanised. After euthanisation, the 
colon tissues of the mice were harvested. Tissues were washed with sterile 1 X 
phosphate-buffered saline to collect faecal contents. The colon tissues were fixed in 
4% paraformaldehyde for pathological evaluation and immunohistochemistry (IHC). 
All animal experiments were approved by the Committee for Laboratory Animal 
Management of Tsinghua University and conducted in accordance with the National 
Institutes of Health guidelines on the ethical use of animals. 

Immunohistochemistry and analysis. IHC was performed using standard protocols. 
The following antibodies were used: anit-p65 (sc-372, Rabbit, Santa Cruz), 
anti-PPAR (ZS-72730, Mouse, ZSGB-BIO), anti-p53 (sc-6243, Rabbit, Santa Cruz), 
anti-COX-2 (sc-376861, Mouse, Santa Cruz). anti-|3-catenin (sc-7894, Rabbit, Santa 
Cruz), anti-CCR2 (ARP 58409-P050, Rabbit, AVIVA SYSTEMS BIOLOGY). 

The expression levels of molecules in IHC digital images were quantified by a 
scoring methods'^. Breifly, the digital images were processed by a color de-convolu- 
tion method to separate the DAB and hematoxylin staining. Scores were calculated 
from the intensity histogram profiles of DAB staining. The trans-nuclear rates of p65 
were measured by ImmunoRatio'^^. 

168 rRNA gene library construction and sequencing. Genomic DNA was extracted 
and purified with the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) and 
used for the preparation of 16S rRNA amplicons. PGR primers targeting the variable 
V3 region were used to amplify the 16S rRNA (forward primer: 5'- 
CCTAGGGGAGGGAGCAG-3', reverse primer: 5'-ATTACGGCGGGTGCTGG- 
3'y'. Both primers were added with lUumina paired-end sequencing adapters, and 
the reverse primer was barcoded with a 6-base code (Table SI). Each reaction mixture 
contained 25 |iL 2X PrimerStar HS Premix (Takara, Dalian, China), 0.5 [iL each of 
the forward and reverse primers (10 |iM final concentration), 50 ng template DNA, 
0.2 mM dNTPs and 1 U PrimeSTAR HS DNA Polymerase (Takara, Dalian, China) 
in a final volume of 50 jil. PGR amplification conditions were as follows: denaturation 
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step at 94'"C for 5 minutes, then 25 cycles of 94^C for 45 s, 50^C for 45 s, and 72 "C for 
60 s, and finally 72"C for 5 min. The amplicons were purified with the QIAquick PGR 
Purification Kit (Qiagen, Hilden, Germany). 10 control samples and 11 treated 
samples were randomly assigned to two flow cells of an lUumina MiSeq genome 
analyser with 151-nucleotide paired-end multiplex sequencing. Human genomic 
libraries from our own laboratory were mixed with the libraries as controls. 

Sequence processing. Sequencing reads were processed with the MiSeq SOP pipeline 
of the Mothur software package*^'*^. Read pairs were assembled as contigs, and then 
contigs with ambiguous bases were excluded. The remaining sequences were aligned 
to the V3 variable region of the 16S rRNA sequences in the Silva database*^, and 
poorly aligned sequences were removed. Unique sequences were identified, pre- 
clustered, and screened for chimeras with UGHIME^". Sequences were classified with 
a naive Bayesian classifier, and sequences that did not belong to the bacterial kingdom 
were excluded (with a 0.8 confidence score). Operational taxonomic units (OTU) that 
had an average percent identity of 97% were picked from the remaining sequences. 

Diversity analysis. Before the diversity analysis, the number of sequences in each 
sample was normalized by picking the minimum number of sequences of samples 
(614247 sequences) and randomly selected this number of sequences from each 
sample. OTU-based analyses of the alpha diversity and beta diversity of the microbial 
communities were conducted with Mothur. In alpha diversity analysis, rarefaction 
curves were constructed based on different OTU levels. Ghao 1 and ACE were 
calculated to estimate community richness. The Inverse Simpson index^* and 
Shannon index^^ were calculated to estimate community diversity. The Simpson 
index-based measure and Shannon index-based measure of evenness were calculated 
to estimate community evenness. In beta diversity analysis, the Jaccard index, Yue & 
Glayton theta coefficient^^ and Bray-Gurtis index^^ were calculated to describe the 
dissimilarity between microbial communities. The results were visualized with non- 
metric multidimensional scaling (NMDS). Analysis of molecular variance 
(AMOVA)^^ was used to test if the dissimilarities between sample groups were 
significant. Taxonomy classification of OTU was made by the Bayesian classifier in 
Mothur with the Greengenes database^^ as reference taxonomy. 

Identification of differential OTUs. The 'edgeR' package^^ was employed to evaluate 
the differences in OTU abundance between control mice and all AOM/DSS-treated 
mice as well as the OTU abundance differences between mice treated with different 
AOM/DSS cycles. OTUs with very low abundance, which were defined as less than 
100 reads, or OTUs that were wholly absent in a total number of samples that was 
greater than the number of samples in the smallest group for comparison were 
filtered. For comparison of control and all AOM/DSS treated mice, an exact test^^ was 
employed; for comparison of control group and groups treated with different AOM/ 
DSS cycles, the generalised linear model likelihood ratio test^^ was performed. Both 
tests used a threshold of false discovery rate (FDR) < 0.05. 

OTU profile clustering and network construction. The Short Time-series 
Expression Miner (STEM)^ was used to cluster the OTU relative abundance time 
course data. The OTUs with very low abundance were filtered as above. The relative 
abundance of the remaining OTUs was calculated for the STEM input. The maximum 
number of model profiles was set to 30, and the maximum unit change between time 
points was set to 2. Significant expression profiles were identified with a false 
discovery rate of 0.05. A microbial network was constructed to reveal the patterns of 
similarity/dissimilarity between profiles. Two profiles were connected if the 
Spearman correlation coefficients of their patterns were significant (P-value < 0.05). 

To construct molecule-OTU dynamic networks, the Pearson correlation coeffi- 
cients between the expression levels scores (trans-nuclear rate for p65) of molecular 
makers and OTU relative abundances were calculated. The network were constructed 
by keeping OTUs which at least significantly correlated with one of the six molecular 
markers (P-value < 0.05). The nodes in the network were coloured according their 
fold changes comparing to the control group. 
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